Geometry dependence of the clogging transition in tilted hoppers 
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We report the effect of system geometry on the clogging of granular material flowing out of flat- 
bottomed hoppers with variable aperture size D. For such systems, there exists a critical aperture 
size D c at which there is a divergence in the time for a flow to clog. To better understand the origins 
of D c , we perturb the system by tilting the hopper an angle 6 and mapping out a clogging phase 
diagram as a function of 6 and D. The clogging transition demarcates the boundary between the 
freely-flowing (large D, small 9) and clogging (small D, large 6*) regimes. We investigate how the 
system geometry affects D c by mapping out this phase diagram for hoppers with either a circular 
hole or a rectangular narrow slit. Additionally, we vary the grain shape, investigating smooth 
spheres (glass beads), compact angular grains (beach sand), disk-like grains (lentils), and rod-like 
grains (rice). We find that the value of D c grows with increasing 8, diverging at n — 8 r where 6 r is 
the angle of repose. For circular apertures, the shape of the clogging transition is the same for all 
grain types. However, this is not the case for the narrow slit apertures, where the rate of growth of 
the critical hole size with tilt angle depends on the material. 

PACS numbers: 45.70.-n, 45.70.Ht, 45.70.Mg 



I. INTRODUCTION 

Granular media are non-equilibrium, athermal systems 
that continue to evade a comprehensive physical descrip- 
tion [IH3] . One of the more intriguing properties of gran- 
ular flow is the phenomenon of clogging, where the flow 
is spontaneously arrested due to obstacles or boundaries 
[4TH"3"] . This can be seen in everyday life, whether as salt 
clogging in a shaker, or pedestrian traffic at a theater 
exit [H]. Understanding and controlling clogging is im- 
portant for industry and agricultural processing. In most 
applications, the goal is typically either to eliminate clog- 
ging from a system entirely, or else to use clogging to con- 
trol output volume, such as metering a bit of salt from a 
shaker. Furthermore, clogging, together with fundamen- 
tal questions about the nature of the glass transition, has 
been responsible for motivating much of the rich held of 
jamming (for reviews, see (151 I16| h 

While most work on jamming has focused on spatially 
uniform systems, clogging behavior is determined by the 
convergent flow toward an opening in the boundary. The 
sudden formation of a stable arch over the exit halts all 
flow, and we can then consider the entire system to be 
jammed. Other systems which demonstrate clogging in- 
clude the flow of vortices through an array of pinning 
sites in type-II superconductors as well as grains flow- 
ing through an array of obstacles [PTl - fH)] . These systems 
demonstrate a clogging transition which depends on the 
density of pinning sites and the packing fraction of par- 
ticles. 

As a canonical example of clogging, we consider the 
discharge of grains from the circular hole of diameter D 
located on the floor of a flat-bottomed hopper or silo. For 
three-dimensional hoppers, it has been observed that the 
time scale one must wait for grains to discharge before a 
clog occurs diverges at a critical aperture size D c [20] 121] . 
We understand this as a clogging transition. For large 



apertures, D > D c , the granular flow will never clog, 
while for the small apertures, D < D c , the flow proceeds 
steadily for a while but eventually clogs. The location 
of the clogging transition is therefore found at D = D c . 
Zuriguel et al. [20] found the value of D c by measuring the 
average mass of grains discharged before a clog occurs, 
(to), as a function of the exit hole diameter D. The value 
of D at which (m) diverges is therefore the critical hole 
diameter D c . They found for spherical grains that D c 
scales with the grain diameter d, as D c /d = 4.94 ± 0.03. 
In their work, they did not see a dependence of D c /d with 
either grain surface roughness or polydispersity. How- 
ever, they did observe that D c /d depends on the ma- 
terial shape. Additional work found that vibrating the 
entire hopper reduced D c [21] . These studies on the clog- 
ging transition for three-dimensional hoppers stand in 
contrast to similar work with two-dimensional hoppers, 
where it was unclear whether the time to wait until a 
clog occurs diverged at a finite hole size [22T 25 . For 
such geometries, the clogging behavior can be equally 
well described by a diverging time scale [221 [23] and an 
exponential [2H [3S] or exponential squared [221 [23] func- 
tions of D. 

The formation of arches over a hole is affected by not 
only grain and hole sizes and shapes, but also by the angle 
8 at which the plane of the hole is tilted away from hor- 
izontal. This was explored by Sheldon and Durian |26j . 
who measured both the flux of discharging grains and the 
clogging behavior as a function of hole size and tilt an- 
gle. In analogy to the start and stop heights for granular 
flow down inclined planes [2 7) , they mapped out a "clog- 
ging phase diagram" by continuously tilting the hopper 
and measuring the angles at which flow spontaneously 
stops (for increasing tilt angle) and starts (for decreas- 
ing tilt angle). Since this gave two critical angles for a 
given hole size, and since both could depend on tilt rate, 
the clogging transition was not precisely located. Fur- 
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thermore, the connection between these start and stop 
angles and the critical aperture sizes of Refs. [20Tf23] was 
not explored. 

In this paper we map out the clogging phase diagram, 
determining a well-defined clogging transition in D — 9 
parameter space for various hopper and grain geometries. 
We do this by observing the values of D where there is 
a divergence in the mass of material discharged between 
clogging events. This is similar to the method used in 
prior experiments, which were only performed for circu- 
lar holes and at 9 = [20 1 [21 ] . Th e value of D c found in 
those experiments describes the location of the transition 
at 6 = 0. We find D c (9 = 0) to be in agreement with 
those values. We measure the variation of D c with 9, and 
discover that the location of this transition is in agree- 
ment with the preliminary measure given by the stop 
angles of Ref . [26] . We also map out the clogging transi- 
tion for a variety of material shapes and sizes (see Table 
|l|), and two different aperture shapes: a circular hole and 
a rectangular narrow slit. Remarkably, for circular holes, 
we find that the shape of the clogging phase diagram is 
identical for all materials. However, this is not the case 
for clogging at a narrow slit. 



II. EXPERIMENTAL DETAILS 

We study clogging for a variety of granular materi- 
als: Potters 2 mm and A100 series glass spheres, U. S. 
Silica Ottawa sand, lentils, and Basmati rice. Table [I] 
lists the relevant properties of these grains, including the 
draining angle of repose 9 r , the bulk density of the gran- 
ular packing p b , and the long d and short d s h or t axes 
of the grains. The maximum tilt angle, beyond which 
the medium loses contact with the wall and no grains 
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can exit, is 9 r , 



We measure 9 r by pouring 



the granular materials into a cylindrical container with a 
large circular hole in the center of the flat bottom. After 
the discharge comes to completion, 9 r is measured as the 
angle which the surface of the remaining grains makes 
with the horizontal. 

To further characterize the grains, we also measure the 
mass discharge rate W versus hole diameter D. The re- 
sults are well-described by the usual Beverloo equation, 

W = C Pb ^{D-kdf/\ (1) 

= C[p b ^g~d?]{D/d-kfl 2 , (2) 

where C and k are dimensionless fitting parameters, and 
g = 9.8 m/s 2 [2SJHH]. For the d = 2 mm glass spheres, 
discharge rate data are divided by p b \/ gd 5 and plotted 
versus D/d'm Fig.[l] along with prior data for two other 
sphere sizes [26]. All three data sets then collapse, as ex- 
pected from the Beverloo equation. The simultaneous fit 
to this equation is very good, as seen in the log-log plot of 
Fig. [iji as well as in the linear plot of Fig. [TJa, where the 
dimensionless discharge rate is raised to the 2/5 power 
thus making the Beverloo form a straight line that van- 
ishes at D/d = k. The fitted value of k for the glass 
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FIG. 1: (Color online) (a) Dimensionless mass flux 
W I (pb\/gd 5 ), and (b) the same quantity raised to the 2/5 
power, for glass beads as a function of dimensionless hole di- 
ameter D/d, where W is the mass discharge rate, pb is the 
bulk density of the packing, g = 9.8 m/s 2 , and d is the grain 
diameter as labeled. The data for the d = 2 mm spheres are 
new, while the other data are taken from Ref. [26]. The fit 
to the Beverloo form, Eqs. l ] 1 12[ ). gives C = 0.56 and k — 1.5, 
and is shown a solid gray curve (a) and line (b). The shaded 
region to the left of D/d — 4.7 indicates where the system is 
susceptible to clogging 



spheres is 1.5 ± 0.1, which is about three times smaller 
than the dimensionless critical hole size D c /d = 4.7 at 
which we find in later sections for the location of the 
clogging transition. Note in Fig. [I] that the discharge 
data is smooth and continuous, in accord with the Bev- 
erloo equation on both sides of the clogging transition. 
This same behavior is also observed for the other grain 
types, for which the fitted values of k are collected in 
Table H 

We investigate clogging through two aperture shapes: 
circular holes and narrow rectangular slits. For the circu- 
lar holes we mount camera irises at the bottom and sides 
of a hopper so that we may explore a continuous range of 
hole diameters D. This flat-bottomed aluminum hopper, 
or silo, has a square cross-section of inner dimensions 
9.5 x 9.5 cm 2 and height 91 cm. For values of 9 less 
than about 60° , the hole at the bottom of the hopper is 
open. For larger 9, we instead use an iris mounted on the 
sidewall at 5.0 cm from the bottom. These two configu- 
rations are equivalent: it has been shown that it is the 
orientation of the aperture, not the container, which in- 
fluences granular discharge [53] . This can also be seen by 
the absence of any discontinuity or kink in the clogging 
transition curve at 9 = 60° (shown later, in Fig. [8]). 

We also investigate clogging for narrow rectangular 
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TABLE I: Material properties and clogging results for the grains tested. The long and short axes are denoted d and d s hort, 
respectively. The bulk density is pb and the draining angle of repose is 9 r . The value k is the dimensionless hole size where the 
Beverloo equation (Eq. predicts the flux to vanish. The dimensionless critical aperture sizes at zero tilt, D co /d, for both 
the circular hole and rectangular slit are also shown below. Note that for all materials D co /d > k. Error bars for d and d 3 hort 
indicate the standard deviations of the distributions. For the other material properties, error bars indicate the standard error 
in the measurements. 



Material 


d (mm) 


d s hort (mm) 


Pb (g/cm 3 ) 




k 


Hole D co /d 


Slit D co /d 


Glass Spheres 


2.02 ±0.04 




1.62 ±0.01 


23.5° 


1.5 ±0.1 


4.7 




Glass Spheres 


0.96 ±0.05 




1.56 ±0.02 


20.0° 






1.7 


Ottawa Sand 


0.77 ±0.10 




1.61 ±0.06 


33.6° 


1.7 ±0.6 


6.0 


2.1 


Basmati Rice 


7.4 ±0.5 


1.48 ±0.10 


0.81 ±0.04 


34.2° 


0.8 ±0.2 


2.6 


1.3 


Lentils 


5.7 ±0.5 


2.3 ± 0.2 


0.83 ±0.02 


30.2° 


1.0 ±0.4 


2.8 


1.4 



slits of constant length 149 mm. For this, a custom-made 
slit of adjustable width D is mounted on the bottom or 
side of a flat-bottomed container. We bevelled the edges 
on the outside of the slit to ensure that the aperture wall 
thickness does not affect the clogging behavior. For 9 less 
than about 60° , we use a hopper with inner cross-section 
28 x 20 cm 2 , height 23 cm, and the slit on the hopper 
floor. For larger values of 9, we use a hopper with in- 
ner cross-section 22 x 14 cm 2 , height 30 cm, and slit on 
the hopper sidewalk The clogging behavior of all four 
classes of materials was tested for both the hole and the 
slit. However, the 2 mm glass spheres were only used 
with the circular hole and the A100 glass spheres were 
only used with the slit. For the 2 mm spheres, we refilled 
the hopper during discharge. For the other materials, we 
did not refill the hopper during flow, and restricted our 
maximum time window to the time to nearly exhaust the 
entire hopper. In all cases, we ensured that the height 
of the grains at the end of discharge was larger than ap- 
proximately twenty grain diameters from the aperture. 
In such a limit, the discharge behavior is independent of 
the filling height. Experiments [35J and simulations [3D] 
confirm that flux is generally independent of filling height 
as long as the height is larger than the hole size. 



III. DISCHARGE DISTRIBUTIONS 

As a prelude to measuring the clogging transition by 
the method of Refs. [2"0T - f2"3"] . now as a function of tilt 
angle, we must first characterize the distribution of the 
mass discharged for a given hopper and grain geometry. 
To initiate flow from a clogged configuration, we break 
the arch over the exit by gently poking it either with a 
wire (for holes) or stiff paper (for slits). By repeating ten 
or more times, we acquire statistics for the cumulative 
distribution, C(m), which is the fraction of events with 
mass less than or equal to m. Fig. [2] illustrates typical 
behavior for glass spheres discharged at zero tilt angles 
through (a) holes and (b) slits of different sizes. After 
scaling by the average discharge mass (m), the data for 
holes collapse nicely to a rising exponential G(vnj (m)) = 
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FIG. 2: (Color online) Cumulative distribution function C 
of mass discharged m, normalized by the average mass (m) 
for a variety of hole diameters and slit widths D as labeled. 
Shown is C(m/(m)) for (a) d = 2 mm diameter glass spheres 
discharged from a circular hole and (b) d = 1 mm diameter 
glass spheres from a rectangular slit, at zero tilt. The overlaid 
curve is the cumulative distribution function for an exponen- 
tial distribution. The distribution for the circular holes is 
similar to an exponential; for the rectangular slit it is some- 
what sharper than an exponential. These distributions are 
typical for all angles and materials studied. 



1 — exp(— m/(m)), exhibiting no trend in hole diameter 
D. The discharge distribution, which equals dC(m)/dm, 
is therefore exponential and it follows that the clogging 
is a random Poisson process that occurs with some rate, 
independent of prior history. This agrees with previous 
observations [2UH25] [SU [32]. For slits, in Fig. [2]b, the 
data do not collapse as nicely for holes but yet display no 
particular trend with slit width. The rise from zero to one 
is somewhat faster than exponential. This differs from 
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FIG. 3: (Color online) Standard deviation a of the discharged 
mass as a function of the average mass (m) for the (a) circular 
hole and (b) rectangular slit. Symbol shapes denote different 
grain types, with shading indicating the tilt angle 9/9 max . All 
aperture sizes D are represented on the plots: for any given 
grain type and 8/6 max , (m) increases with hole diameter or 
slit width D. The standard deviation is roughly proportional 
to the average. The expected relationship for an exponential 
distribution is a — (m), overlaid. The discharge distributions 
through a slit are somewhat sharper than expected for an 
exponential distribution. 



TD 

< 

Q. 





a 


hole 




10 4 






// /' - 


10 3 








10 




r' // 


r 

^y' o Glass spheres 


10 






-Jffl a Sand 




a ,4 y^ 


®-J^ □ Rice 
x Lentils 


1 






— <m> = F/(D C -D) Y 


0.1 




■' I 


- - - <m> = F Exp(B D 2 ) 
i i i 



1 2 3 4 5 



TD 
< 

7\ 

E 
v 



1000 


I I 
b) slit 


I 


1 1 , / 1 VI: 1 1 /,! 


100 






a 4' 1 4 y f 

/ J ,' /_.' 


10 




ffl 41' 


it %-' 


1 








0.1 


i 


1 


1 1 1 1 1 



0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 

D/d 



FIG. 4: (Color online) Dimensionless measure of the average 
mass discharged (m) / (pt,Ad) before a clog occurs vs. dimen- 
sionless aperture size D/d. A is the aperture area, pi, is the 
material bulk density, and d is the long grain axis. For the 
circular hole, D is the hole diameter, while for the rectangu- 
lar slit, D is the width of the slit. Plotted are data for (a) 
the circular hole and (b) the rectangular slit, for 6 = 0. The 
average mass grows with D faster than exponential squared 
(indicated by the dashed curve). Instead, (m) vs. D is best 
described by the diverging power law of Eq. |3| , indicated by 
the solid curves. Eq. (J3J diverges at the critical aperture size 
D c , indicated by the vertical dotted lines. 



a previous experiment for a wedge hopper with a long 
narrow slit [33], where flow was initiated by jolting the 
hopper and where a power-law distribution was observed. 

The standard deviation a of the discharge masses is 
a simple parameter for characterizing the width of the 
distribution. Results for a for all grain types and tilt 
angles are collected in Fig. [3] for (a) holes and (b) slits. 
For holes, the data span many orders of magnitude and 
closely agree with a — (m), which is the expectation 
for a Poisson process with an exponential distribution. 
For slits, a approach (m) for large (m), i.e. for large slit 
widths. For smaller slit widths, a < (m) holds and hence 
the discharge distribution is sharper than exponential, as 
seen explicitly for the glass beads at 8 = in Fig. [2Jd. In 
all cases, we thus find that the distribution of discharge 
masses is exponential or sharper. Therefore the average 
mass (to) can be determined to within ±30% uncertainty 
from only ten trials and used, next, to locate the clogging 
transition. 



IV. CLOGGING TRANSITION 

We now determine the location of the clogging transi- 
tion by using the measurements of the average discharge 
mass (to) for different values of the opening size D. Fig. [4] 
shows how the discharge mass grows with aperture size 
for the various materials when 9 = 0. For comparison 
between the different materials, we have scaled (to) by 
PbAd, where A is the aperture area, pb is the material 
bulk density, and d is the long axis of the grains (see Ta- 
ble [I]) . As such, the y-axis represents, in units of d, the 
typical height of a column of grains over the aperture 
that is discharged before a clog forms. As observed in 
the figure, the average discharge mass (to) grows with 
D, and does so faster than exponential squared (indi- 
cated by the dashed curve). In fact, the behavior can 
be well described as a power law of D c — D, with (to) 
diverging at a critical aperture size D c : 

(m) = {Dc F _ D) -n (3) 
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FIG. 5: (Color online). Fitted exponent 7 from Eq. Q 
vs. 6/6 max , when 7 is allowed to float as a free parameter 
in the least-squares fitting of average mass discharged (m) 
vs. aperture size D. The error bars 5j indicate the fitting 
confidence intervals for 7. The average value of 7, weighted 
by l/(Sj) 2 , is constant for values of BjQmax < 0.8. These 
values are 7 = 5.1±0.3 for the circular hole and 7 = 1.9±0.1 
for the slit. Therefore, the values of 7 = 5 and 7 = 2 are used 
as fixed parameters when fitting Eq. Q to the data, indicated 
by the solid lines. Some of these fits are shown in Fig. [6] 



FIG. 6: (Color online) Average mass discharged vs. aperture 
size D for glass spheres clogging in the (a) circular hole and 
(b) rectangular slit. The average mass increases with aperture 
size, diverging at a particular critical aperture size D c . Fits 
to the diverging power law of Eq. Q are overlaid, with the 
exponent 7 = 5 for the circular hole and 2 for the rectangular 
slit, respectively. This behavior is typical for all angles and 
materials studied. 



where F, D c , and 7 are fitting parameters J20-23J. This 
contrasts with two-dimensional hopper discharge, where 
the average discharge mass or duration can be fit well by 
Eq. (jsj, but can equally well be described as growing with 
D exponentially [2H L2S] or as an exponential squared 
[22), [23] . These forms fail to be generally consistent with 
all the materials shown in Fig. [4j although all forms fit 
the rice data. 

For our system, we can determine the value of D c by 
fitting the data to Eq. These fits are overlaid as the 
solid curves in Fig. [4] With three free fitting parameters, 
the fit does not converge for the rice clogging at 9 = 0. 
This difficulty occasionally occurs for other materials at 
various values of 9 as well. Therefore, we next describe 
our method for determining the best value of 7 to use for 
performing these fits and finding D c {9). 

To determine the values of 7 to use in fitting, we first 
fit the data for each material and 9 to Eq. ^ , allowing 7 
to float as a free fitting parameter. Those fitted values of 
7 for which the fits converge are shown in Fig. [5] The er- 
ror bars Sj indicate the fitting confidence intervals. The 
exponents are scattered over a range of ±2 for the hole 
and ±4 for the slit. To within error bars, 7 appears con- 
stant for 9/9 ma x < 0.8. The weighted average value of 



7 over that range is 5.1 ± 0.3 for the hole and 1.9 ± 0.1 
for the slit. The weighted average for all 9 is 3.8 ± 0.2 
for holes and 1.4 ± 0.1 for slits. Based on these results, 
we now repeat the fits using fixed values of 7 = 5 and 
7 = 2 for the holes and slits, respectively. This gives 
good fits to all the materials and values of 9 for which 
we have data. For example, Fig. [6] shows (m) vs. D for 
selected values of 9 for the glass spheres, overlaid with 
fits to Eq. ([3]). This plot is typical for all the angles and 
materials studied. 

As a final remark, we note that Zuriguel et al. found 
7 = 6.9 ± 0.2 for a three-dimensional hopper [20]. For 
a two-dimensional hopper, it has been observed that 
(to) vs. D could be effectively described by both an 
exponential-squared function of D and by Eq. ^ [2"2"ll23| . 
Fits to data for two-dimensional hoppers have produced 
values of 7 = 11.2 [22 j and 7 = 12.7 ± 0.1 [23]. Ad- 
ditionally, Mankoc et al. found values of 7 = 7.6 ± 0.5 
and 7 = 8.6 ± 0.2 for vibrated and non-vibrated three- 
dimensional hoppers, respectively [21) . Fortunately, the 
actual choice of 7 is not critical for our purposes: chang- 
ing 7 over the range of 1.5 — 9 does not significantly 
change the D c values to be presented in the next section. 
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FIG. 7: (Color online) Clogging phase diagram for the glass 
spheres for the (a) circular hole and (b) rectangular slit. The 
squares, shaded by the average discharge mass (m), indicate 
locations in parameter space where data was collected. The 
solid circles show the location of the critical hole diameters 
or slit widths D c , as determined in Fig. [6] These values of 
D c describe a well-defined transition between the clogging and 
freely-flowing regimes. As indicated by the dotted and dashed 
lines, respectively, the transition for the circular hole is near 
the angles where the flow has been observed to spontaneously 
start and stop during continuous tilting of the hopper [26] . 
Fig. [8] displays the clogging transitions for all the materials. 




FIG. 8: (Color online) Clogging transition curves for all ma- 
terials tested. We scale the tilt angles 8 by the maximum 
possible tilt angle for the materials f? max = w — 8 r . We also 
normalize the critical hole diameters and slit widths D c by 
dividing by D co , the value of D c at 8/8 max = (see Table [I] 
for values of 8 r and D co /d for the various materials). The 
transition occurs at the same location in parameter space for 
all materials flowing through the (a) circular hole. However 
for the (b) narrow slit the transition depends on the material. 
The relations given by equations < |4|5[ ) are overlaid. Both rela- 
tionships describe the clogging transition for the circular hole 
well, but not the slit. 



V. CLOGGING PHASE DIAGRAM 

We now map out clogging phase diagrams, which spec- 
ify whether a given material flows freely or is susceptible 
to clogging, as a function of aperture size and tilt an- 
gle. Specifically, we follow the previous section in using 
Eq. ([3| to fit for the critical hole sizes D c at which the 
discharge mass diverges for a set of different tilt angles 
ranging from zero to # max = n — 9 r . For glass spheres, the 
results are collected in Fig. [7] for both (a) holes and (b) 
slits. This figure displays a field of data points where, for 
a given {9, D}, the value of the average discharge mass 
is indicated by the shading. For a given tilt angle, the 
color darkens as D increases and the divergence at D c is 
indicated by a solid circle. The locus of D c values, for 
different tilt angles, are joined by solid line segments and 
thus serve to separate the 6 — D parameter space into 
two regions where clogging does and does not occur. 

The qualitative shape of the clogging phase diagrams 
in Fig. [7] is the same for both circular holes and narrow 
slits. For both, the transition rises steeply since D c is 
found to be nearly constant for small tilt angles. For 
large tilt angles, greater than about about ninety de- 
grees, D c increases rapidly and, intuitively, diverges at 



$max = it — 8 r . The shape of the transition for cir- 
cular holes may be compared with the earlier estimate 
from Ref. [26] based on start and stop angles. As seen 
in Fig. [7^,, the locus of critical hole sizes matches quite 
nicely with the locus of stop angles. While this reinforces 
the validity of Ref. [26J , more importantly it shows that 
the clogging phase diagram for a specific material and 
aperture geometry may be confidently characterized in 
terms of stop angles, which are far easier and faster to 
measure than critical hole sizes. 

Next we investigate how clogging phase diagrams such 
as Fig. [7] are affected by the geometry of the grains. As 
above, the locus of critical hole sizes D c are found at a 
set of tilt angles, now for sand, rice, and lentils. The vari- 
ous grain types have different sizes and different angles of 
repose, which affect the location of the locus of clogging 
transitions in the 9 — D plane. Therefore, to scale this 
out and facilitate comparison, we normalize the aperture 
size by D co , the critical hole diameter or slit width at 
9 = 0, and we normalize the tilt angle by max = tt — 6 r . 
Values of D co are given in Table [TJ The resulting scaled 
clogging phase diagrams for the four grain types are col- 
lected in Fig. [8] for (a) holes and (b) slits. Remarkably, 
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we find a very good collapse of the transition data for 
all four grain types in Fig. [8^, for circular holes. This 
suggests that the shape of the clogging phase diagram is 
universal, independent of grain type. But by contrast, 
for long narrow slits, the transition curve data in Fig. [SJa 
show no such collapse for the different grain types. The 
universality is only for circular holes. 

At present we know of no theory that can predict the 
universal shape of the clogging transition for circular 
holes. We believe that this is an important challenge, 
in the same spirit as the notion of deep commonality 
for jamming transitions. For clogging, this commonality 
now takes quantitative form. To facilitate comparison 
with future theories, we note that the clogging transition 
in Fig. [8^, can be very well described empirically by the 
equation 

where we find j3 = 0.56±0.02. This function is plotted as 
a black solid curve in Fig. [8| For (a) holes, it matches all 
the data quite well. For (bj slits, it matches the data for 
lentils but is quite wrong for the other grain types. For 
small angles (D/D co ~ 1), Eq. Q may be expanded and 
solved for 9/0 max = y/(D/D~- 1)/(1 -/?). Another 
potential form of clogging transition for circular holes is: 

C0S (^7T^) =( D c»/D) 2 . (5) 

\ ^ "max J 

This equation is also plotted on Fig. [8] as a grey dashed 
curve. The behavior of this form is very similar to Eq. Q , 
especially at small 9. Here, Eq. ^ can be expanded for 
small 9 (D/D co ~ 1) as 0/9 max = (2/n)^D/D co - 1. 

VI. CONCLUSION 

In this work, we have reported a well-described clog- 
ging transition for the discharge of noncohesive grains 
from a hopper, which we find by observing how the av- 
erage mass discharged between clogs, (m), grows with 
aperture size D. Consistent with earlier work, (m) grows 
as a power-law of 1/(D C — D), with an exponent 7 that 
depends on the aperture shape: 7 = 5.1 ± 0.3 for cir- 
cular holes and 7 = 1.9 ± 0.1 for slits. This indicates 
that the hopper exit geometry is key in determining the 
clogging behavior. These values are somewhat smaller 
than those reported in the literature for two-dimensional 
hoppers [?:]] and three-dimensional hoppers with cir- 
cular holes [201 HI]- However, we are aware neither of 
prior work reporting 7 for slits, nor of any models for 
predicting the values of 7. 

Using the critical aperture sizes D c , we map out a 
clogging phase diagram as a function of hopper tilt an- 
gle 9 and aperture size D. We have shown that the 



phase diagram can be scaled by the two system param- 
eters m ax and D co . These parameters, together with 
the shape of the clogging transition curve, completely de- 
scribe the clogging phase diagram. The clogging behavior 
at 9 = agrees with similar experiments , exhibit- 

ing a well-definied clogging transition with finite critical 
aperture size, which defines D co . However, we now have a 
clear picture of the shape of the clogging transition curve 
as the hopper is tilted and 9 grows. This result also is 
in agreement with the clogging transition sketched out 
by Ref. for glass spheres and a circular hole. And, 
significantly, the shape of the transition curve is univer- 
sal for circular holes, but depends on the nature of the 
grains for long narrow slits. 

Besides elucidating the value of 7 and the universal 
shape of the clogging phase diagram for circular holes, an- 
other important theoretical challenge raised by our work 
concerns the relation between the Beverloo cutoff hole di- 
ameter, Db ev = kd, and the critical diameter D co for the 
clogging transition at zero tilt angle. As evident from 
Table |lj the ratio D co jD\, ev is noticeably greater than 
one and depends on more than just some power of /i. By 
comparing the glass beads and sand, for which the grain 
aspect ratio is one, we can estimate the \x dependence. 
By scaling this out and plotting vs aspect ratio, we can 
estimate the dependence on long and short grain axes. 
Altogether, this gives 



Dbev/d — 


lV' a (<Wt/d) - 8 , 


(6) 


D co /d = 


7.6»°- 6 (d short /d) 0S , 


(7) 


DcojDbev — 


^■ z {d short /df A . 


(8) 



While these formulae reasonably describe our data, we 
should emphasize that the variation of fi and aspect ratio 
is limited and not independent; furthermore, we have no 
evidence that the intermediate grain axis does not play 
role. Simulations could be very valuable for systemati- 
cally varying grain parameters and establishing the ex- 
tent to which the Beverloo cutoff and the clogging critical 
hole diameters obey some such formulae. Another line a 
future research, which would be particularly important 
for industry, would be to similarly study the clogging of 
grains which are somewhat cohesive. Answering these 
questions regarding the clogging transition will provide 
great insight into understanding the physics of granular 
materials, as well as other far-from-equilibrium and dis- 
ordered systems. 
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